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We reexamine the putative Peierls transition in a (5,5) 
metallic nanotube. We show that the conduction electrons 
at the Fermi level do not couple to the longitudinal acoustic 
phonon but rather to a folded-in graphene zone edge phonon 
having the proper Kekule modulation symmetry. The calcu- 
lation for the mean-field transition temperature gives 15 K, 
a value comparable to previous estimates. We discuss the 
significance of this transition temperature. 

PACS numbers: 71.45.Lr, 63.22.-|-m 

Since their discovery by Ijimatl in 1991, carbon nan- 
otubes have attracted a lot of interest due to their un- 
usual geometry and theiij-structural and electronic PWfj- 
erties. Band calculationsB (confirmed by experimentajc) 
have stressed the one-dimensional (ID) character of 
single- wall carbon nanotubes (SWCNT). They showed 
moreover that, depending on their diameter and he- 
hticity, a SWCNT can either be metallic or semi- 
conducting. The metallic tubes provide an excellent en- 
vironment for one-dimensional (ID) conductor physics. 
The existence of a putative Peierls transition due to the 
electron-phonon coupliiig has already been addressed by 
MintmircQ and HuangQ who predicted the occurrence of 
a Peierls transition at IK and 9K respectively. We pro- 
pose to reexamine the Peierls transition owing to what 
we believe to be a wrong selection of soft phonon mode 
by Huang. The longitudinal polarization of the phonon 
they considered, taken fromO, cannot lead to a Kekule in- 
stability natural to graphene, the two-dimensional (2D) 
carbon sheets that are rolled up to form the tubules. We 
will consider in the following the (5,5) armchair SWCNT 
which presents metallic properties. We shall reevalu- 
ate its phonon spectrum, deduce the electron-phonon 
coupling amplitude, and finally estimate the mean-field 
Peierls transition temperature. 

We first recalculate the phonon spectrum of graphene 
by the dynamical matrix method using the force con- 
stants obtained by MaedaQ. We solve numerically the 
secular equation for the phonons. The polarization vec- 
tor of the phonons is written as 

e{~t,iy) ^ {eA{'q',t^),eB{l',i^)), (1) 

where e,, {if, v) is the polarization of one of the two 
atoms r — A, B in the unit cell for the phonon of wave 
vector (f and branch v. Here e^. has two components along 
the mutually orthogonal directions x and y (the tube axis 
direction) in the sheet respectively (see Fig. 1). These 
polarizations satisfy the relation e {~lf, v) ^ e {if, v) . 
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FIG. 1. Unfolded unit cell of (5,5) nanotube (a) and Bril- 
louin zone of graphene (b). 

There are two phonon groups relevant to our study. 
The first one is at the K points in the 2D hexago- 
nal first Brillouin zone (BZ) of graphene (see Fig. 1). 
Its phonons have the wave vector if — (0, ±47r/3a) 
or (±27r/A/3a, ±27r/3a). The polarizations for if = 
(0, 47r/3a) are given by 
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e(^,l)--(i,l,-i,l),e(^,2) 

e("9',3) = -{-i,-l,~i,l), e{lf,A) = - (i 



(-1,1, i,l), 



4,z,l) (2) 



The second group, to which the acoustic phonon with 
longitudinal polarization belongs, has polarizations 

e(^,l) = -i=(0,l,0,-l),e(^,2)-i=(~z,0,i,0) 
e(^,3) = -^(0,1,0,1), e(Y, 4) = -^ (*,0,i,0) (3) 

for if — (0, — 27r/3a). Here a/\/3 is the carbon-carbon 
distance in graphene. 

The nanotube is constructed by. rolling the sheet of 
graphene in the prescribed fashiona. The (5,5) SWCNT 
has the right diameter to be capped at the ends by half 
a Cgo sphere. Its unit cell has length a and contains 
20 atoms (see Fig. 1). The requirement of periodicity 
around the nanotube imposes the following quantization 
condition on qx = (27r/5\/3a) Tn, where —5 < ttj < 5 is 
an integer. The equilibrium position of an atom on the 
nanotube is specified by {n,s), n indicating the position 
Rn of the unit cell along the nanotube axis and s referring 

to the position dg of the atom within the unit cell, so that 



Rn,s — RnV + ds 



(4) 



Because of the difference in the length of the unit cell 
of SWCNT (5,5) with the one of graphene, the first Bril- 
louin zone has a smaller length in the y direction given by 



It: I a. We can than approximate the phonon spectrum of 
the nanotube by folding the spectrum of graphene, sub- 
ject to the quantization rule for q^. above, into the BZ of 
the nanotube. We can now write the polarization vector 
Ss (g, v^ m) of atom s for the phonon of wave number q 
(along the axis of the tube) corresponding to branch i' 
and quantum number m in terms of those in graphene: 
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Si(2) (q, V, m) = -j=eA[B) i~t, V 



83(4) (g, V, m) 



gi(9+G,)§+ig^ 



where qy = q + Gq for \q\ < ir/a, and Gq = ng(27r/a) is 
a reciprocal lattice wave number of the nanotube. Also, 
one has the property Ss+4 = exp{iqxV^a) Ss- The result 
of this folding- in for the branches of interest here is shown 
in Fig. 2. 

The electron-phonon interaction Hamiltonian, treaied 
within the tight-binding approximation, is written agj 



is now possible to calculate the electron-phonon coupling 
amplitudes for these two conduction bands bw using the 
eigenvectors U{k) given in the article of JishiQ. 

It is important to realize at this point that the momen- 
tum conservation condition q — k — k' — G imposes that 
G 7^ whenever |A;| = |fc'| ~ kp since |fc — fc'| = Air /3a 
is outside the first BZ of the nanotube and must be 
brought back within by a jCj = 27r/a. As a consequence 
of this, one has \q\ = 27r/3a. There are two possibili- 
ties of interest here. The first corresponds to a graphene 
phonon having \qy\ = \q + Gq\ ~ 47r/3a for ?7i = or 
eA(B) {'^-.v) , (5) i^^i ^ 1^1 = 27r/3a for m = 5. These are the K point 
phonons of graphene. The other is \qy\ = \q\ = 27r/3a for 
TO = to which the acoustic phononj-with longitudinal 
polarization used by Huangjj and JishiQ belongs (see Fig. 
2). 



), 



H, 



I— rj Z^ 9k,p-k',p'\^)^k,p,a^k',p', 






(6) 



where L is the length of the tubule, Cjf" ^ {Ck,p,a) is the 
creation (annihilation) operator for an electron with wave 
vector k and spin a in band p, G is a reciprocal lattice 
wave number such that q — k — k' — G is within the na- 
totube's first BZ, and bq^^m is the phonon operator. The 
electron-phonon interaction matrix element g^'™^,/ „i{G), 
is given within the tight-binding approximation by the 
following expression : 



9k',p;k',p'\^r 






hu;k-k'-G,uj 

[VaJiO,s;n\s')]oUlpik)Us',p'ik') 



a,nf,s,s/ 

X {e-^'='«"'[S,(fc-fc'-G,i.,TO)]„ 
-e-*^-^"'[S..(fc-fc'-G,z.,TO)]J, 



(7) 



where /x is the mass per unit length of the tubule, cOq.ij^m is 
the phonon frequency, J{n, s;n , s ) is the matrix element 
of the atomic potential connecting orbitals on nearest- 
neighbor atoms (n, s) and (n', s'), and U{k) is the Bloch 
eigenvector that diagonalizes the purely electronic part of 
Hamiltonian. The elements of U{k) obey to the periodic 
boundary condition Ui.p{k) = f/2i,p(fc) and correspond 
to the graphene Bloch eigenvectors folded into the nan- 
otube's BZ. The Fermi level is situated at kp — (27r/3a) 
at the crossing of two bands having the dispersion rela- 
tions ep(fc) = ±Jo [1 - 2cos(fca/2)] with p == 1,2. The 
band crossing corresponds to the K point of the BZ of 
graphene of Fig. lb. The wave functions at the six K 
points bear the signature of the Kekule bond arrange- 
ment of graphene which can be obtained by coupling to 
the -^ = (0,±47r/3a), (±27r/\/3o, ±27r/3a) phonons. It 
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Wave vector (71/a) 

FIG. 2. Phonon dispersion of the (5,5) nanotube com- 
ing from the folded-ing = {0,qy) and the q = {2Ti/?,^a,qy) 
phonons of graphene which coalesce at g = 27r/3a. The cir- 
cle contains the mode that couples the electrons at the Fermi 
level. The square shows the longitudinal mode previously 
used by other authors. 

It should be apparent by now that the K point phonons 
naturally couple to the Kekule structure whereas the 
acoustic phonons, having strictly longitudinal polariza- 
tion, cannot couple to the transverse (around the nan- 
otube) part of the Kekule modulation. Indeed, the 
circular polarization (order parameter) associated with 
the \qy\ — |g-t-Gg| — At: /3a and m — phonon, as 
can be deduced from the polarizations above, corre- 
sponds to a Kekule modulation. To show this, let us as- 
sume a non-zero order parameter (h^kp-i..- 
bKAie^'^ giving rise to a displacement u„ s — 

^yhLUkF.i,a/'2. Ss(-fc_F,l,0)exp[-i(fcF-Rn - 0)]fe.;fA/ + 
c.c. One obtains the following relative displacements 



dn,3,n,4 = ^D (1, 0) COS [a (n)] , 

dn,2,n,3 = D (cos {tt/3) , siu (7r/3)) cos [a (n) + tt/3] , 
d„,i,„,2 = -D (1, 0) cos (a (n) + 27r/3) , 
d„_i,3,n,2 = D {- cos (7r/3) , sin (7r/3)) cos [a (n) + tt] , (8) 



"kF,l,0 



where D — 2yJhuj-kp.ifl/'2.Q bKM, and a{n) = 
(n — 1) 2tt/3 + ■K 12 — (j). This is the Kekule modula- 
tion with period An = 3 or wavelength Ay — 3a. The 
relative displacements lie along the carbon-carbon bond 
directions leading to maximum modulation of the elec- 
tron hopping amplitude Jo and thus maximum electron- 
phonon coupling. 

In the situation were k ( k') is near +fci? {—kp) is, for 
the K point phonons having tti = 0, 5, one finds 



9k,l;k',2\~~) ^ 9k,2;k' ,1\~~) 



Vl/ = 1, 



(9) 



for p ^ p' and 

,2-K^ 



4,m /^'^s 
■Jk.p\k' .p\ ' 



yk,p;k' ,p\ -^ ^k.p;k' .p 



^k,p;k' ,p^ 
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27r 
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(10) 



for p — p' — 1,2. Here rji, — qoJoy^ha/2mcUJo 
(luo — ^±2-K/3a.i^,o), where rric is the mass of a carbon 
atom, and J{r) is assumed to depend only on distance 
such that qoJo = [9J(|r|)/9 |r|]QQ. It can be seen that 
5fe;^fe',i(27i-/a) = and gl'^]k' .iC^^ / (^) = -6«??i/V40 at 
k = —k' = kp. 

Futhermore it can be demonstrated that 

9lik'./i) -9lik'. 2i^i)-0 V.= l,...,4. (11) 

for all m = and \q\ = 27r/3a phonons (this includes 
the acoustic phonon). This result is different from the 
estimate of the other mentioned authors. We recover 
their result when evaluating g^'^.j./ ^(0) or, equivalently, 
assuming the longitudinal polarization corresponds to a 
K edge phonons. This is probably at the origin of the 
confusion. 

One can now estimate the mean-field transition tem- 
perature. After integrating the phonons in the interac- 
tion representation, the partition function of the system 
can be expressed as 



Z = Z^,Jt L-^""'^ Tr 



./^fH,(r-r')<r-r') 



(12) 



where iJo,e = T.p.k,a ^p (k) Clp^^Ck,p,a, and Hi (t - t') 
is the retarded electron-electron interaction given by 



Hi{t 



E 

p,k,a;p' ,k' ,cr';q 



Do{q,T-T') 



^ 1,0 /27r 
^ 9k,p-k' ,p\~ ) 



1,0 /^^, 

9k,p':k' ,p'\ ) 



^ ^k,p,a^k~q~^,p,a'^k' .p' ,a'^k'+q+^,p' ,g'- (13) 

This last expression can be alternately written as 

Hi (t -t')= Y. 9{r ~ r')OKM{q,T)0*KMiQ,r) 

p.k^{7\p' ,k' ,a' \q 



where 
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g[r-T') = Do{- — ,T-T 



9kF,l--kF,l\ ) 



1,0 



(14) 
(15) 
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KM 



("■-J-T^E 
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k.a 



C+, (k) Ci,, [k-q-'^ 



27r 



-ClAk)C2,a[k-q- — 



(16) 



and the Matsubara time Fourier transform of the phonon 
propagator is DQ{q,ujM) = -2ujq,i^o/ {uj^^i^o + ujIj). The 
order parameter associated with Okm, the operator 
that is naturally selected by the electron-phonon cou- 
pling, is the electronic density Kekule modulation oper- 
ator. To verify this, let us calculate the electron density 
{^*if)^{f)) where ^(r) = Efe,p,.,n,. e''=^"C^a,p(fc)0(^" " 
Rn.s) Ck,p,c!, (t^i^ ~ Rn.s) is the TT Orbital of the car- 
bon atom at Rn,s^ and the average is taken over 
a broken symmetry situation having < Okm >= 



-(cz ,„c_ 



kF,2,ij 



Pkms.^'^ ■ One can show that 

(^*(f)*(f)) = - Y^ K(n,s;n',s') 



Pkm 



X (/)* (f — R„ 



')(t>{r - Rn,s) 



(17) 



in which only the bond densities contribute (n, s ^ n' , s'). 
One has K(n, 3; n, 4) = ^^ cos [n27r/3 — (^] , «;(n, 2; n, 3) = 



-^cos[(n-f l/2)27r/3-^], 
10 ■ 



K(n, 1; n, 2) — ^ cos [in + 1) 27r/3 — Lp] , and «;(n, 2; n — 
1,3) = — :^ cos [(n -I- 3/2) 27r/3 — 93] . Comparing these 
bond density modulations to the relative phonon dis- 
placements calculated earlier, one sees that the two co- 
incide when putting cj) = ip + it/2 — 27r/3. 

We can now calculate the static Kekule modulation 
susceptibility for the system at q = 27r/3a. The random- 
phase approximation (RPA) result is given by : 



Xkai 



X 



l~g{ujM^O){xl + xl) 



(18) 



where Xo iXo) is the bare RPA charge response func- 



tion of the first (second) band , xk 



M 



X 



X , and 



X = Xo + Xo - 25(0)xJXo -As 5(0) and Xo + Xo are 
both negative, the susceptibility can diverge when the 
denominator goes to zero. This correspond to a mean- 
field critical temperature given by 



Tc = 3.86 Joexp 



1 



2N{eF)g{Q) 



(19) 



To obtain this pre-exponential constant, we did a nu- 
merical calculation of Xo (Xo) using the exact dispersion 
given above and fitted it to a standard Cln(D/T) form. 
Here, one has N (ep) = 2/ (n^/SaJo) and ^(0) = — y^ 
((7q Jg a) / (mcCjQ^). For a numerical estimate of Tc, 

we take! a = 2.46 A, Jq = 2.2 eV, Qq = 2.5 A~\ along 
with our value of Hujq = 0.166 eV corresponding to the 
graphene K point phonon. We obtain Tc — 15.2 K. This 
estimated Peierls temperature is of the same order as 
Huang's even though the selected phonon mode is dif- 
ferent. The reason is that although the numerical factor 
for our effective electron-phonon coupling amplitude is 
larger, our phonon frequency is also larger. The two ef- 
fects nearly cancel each other. 

Even though nothing much has changed for the nu- 
merical estimate of the Peierls-Kekule mean- field transi- 
tion, we feel it is important that people use the correct 
electron-phonon interaction coupling to the right phonon. 
This is important in conductivity or superconductivity 
calculations. 

The predicted Tq would of course change if one were to 
take into account the modification in stiffness constants 
due to the curvature of the SWCNT. These changes 
would be quantitative, not qualitative. One may of 
course, because of the Mermin- Wagner theorem, ques- 
tion the validity of a mean- field transition in ID sys- 
tems. There is of course no transition. Ttue mean-field 
Tc signals the occurrence of a pseudo-gapO. Moreover, 
the mean-field formula is for adiabatic phonons. This is 
surely not the case here as the thermal energy at the pu- 
tative transition kgTc is much smaller that the phonon 
energy. Quantum fluctuation effects would be consider- 
able. They would wipe out the pseudo-gap. 

The above analysis also presumes of an infinitely long 
nanotube. These are of course of finite (micron scale) 
length. The analysis would prevail as long as the coher- 
ence distance is smaller than the tube's length. 
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